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INTRODUCTION 


Solid  rocket  motors  are  considered  stable  if  the  acoustic  losses  within  their  interiors 
are  larger  than  the  acoustic  gains.  The  major  sources  of  acoustic  gain  and  loss  are, 
respectively,  the  unsteady  combustion  process  and  the  convection  of  acoustic  energy 
through  the  nozzle.  In  a  series  of  pioneering  analytical  studies,  Culick14  identified  the 
interaction  between  the  transverse  injection  of  mean  combustion  gases  into  the  rocket 
interior  and  the  longitudinal  standing  wave  sound  field  as  a  secondary  but  nonetheless 
important  source  of  acoustic  loss.  Because  the  interaction  takes  place  within  the  motor 
interior  where  the  injected  flow  undergoes  continuous  turning  prior  to  its  convection  out  the 
nozzle,  these  losses  have  been  identified  as  acoustic  "flow-turning"  losses.  Culick 
observes  that  during  combustion,  exhaust  gases  are  injected  into  the  rocket  interior  with 
initially  zero  axial  acoustic  momentum  and  exit  through  the  nozzle  with  nonzero  axial 
acoustic  momentum.  He  argues  heuristically  that  this  increase  in  momentum  must  be 
acquired  at  the  expense  of  an  irreversible  loss  of  acoustic  energy  by  the  mean  gases  within 
the  rocket.  Despite  its  widespread  use,  Culick's  theoretical  model  appears  to  be 
misunderstood  by  many  members  of  the  combustion  community. 

The  objective  of  this  investigation  is  to  examine  Culick's  model  both  experimentally 
and  theoretically.  The  intent  here  is  to  clarify  apparent  misconceptions,  based  on  a  survey 
of  the  open  literature,  concerning  its  interpretations  and  limitations.  A  survey  of  the 
literature  also  revealed  the  need  for  further  data  to  support  Culick's  claim  that  the  injection 
of  fluid  into  a  chamber  containing  longitudinal  sound  results  in  a  net  absorption  of  sound 
and  is  hence  stabilizing. 

In  Reference  1,  Culick  investigated  the  stability  associated  with  the  injection  of  fluid 
into  a  chamber  containing  longitudinal  sound.  Although  the  model  was  one-dimensional, 
his  investigation  was  quite  general  in  that  it  included  the  effects  of  combustion,  particulate 
matter  and  variable  chamber  area.  Since  the  interest  here  is  on  understanding  the  basis 
mechanics  identified  with  "flow-turning"  losses,  these  effects  will  not  be  investigated  in 
this  report.  Instead,  the  far  simpler  problem  is  investigated  of  understanding  "flow¬ 
turning"  losses  for  the  case  of  injection  of  ambient  (cold)  air  into  a  channel  containing 
longitudinal  sound.  Culick  investigated  this  problem  and  derived  the  "flow-turning" 
expression  shown  below  using  the  following  assumptions: 

(i)  the  gases  are  inviscid,  non-heating  conducting  and  obey  the  perfect  gas  law; 

(ii)  the  gases  were  transversely  injected  into  the  chamber  at  very  low  speeds; 

(iii)  the  chamber  was  closed  with  rigid  end  walls; 

(iv)  the  chamber  contained  longitudinal  sound. 

For  the  interested  reader,  a  detailed  derivation  of  Culick's  "flow-turning"  losses  is 
presented  in  the  Appendix.  Culick4  derived  the  following  expression  for  the  decay 
constant  (otpy)  due  to  the  "flow-turning"  losses. 
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Equation  ( 1 )  shows  that  the  interaction  between  the  transversely  injected  mean  flow 
(mb/P  =  v  b )  —  which  is  modelled  as  a  source  of  mass  at  the  chamber  boundary  in  a 
one-dimensional  formulation  -  and  the  longitudinal  standing  wave  pattern  within  the 
chamber  results  in  a  loss  of  acoustic  energy.  Culick2  claims  that  the  "flow-turning" 
expressions  (apj-)  "are  approximate  representations  for  the  influence  of  essentially  viscous 
processes  in  the  real  flow.  The  viscous  processes  occur  in  thin  regions  near  the  surface 
and  are  present  when  there  is  a  component  of  unsteady  motion  parallel  to  a  surface  through 
which  mass  is  entering  or  leaving."  Culick  apparently  believes  that  all  the  "flow-turning" 
losses  take  place  within  this  region.  Although  Culick  never  explicitly  defines  this  "thin 
region",  it  is  certainly  reasonable  to  presume  that  it  is  related  to  the  local  acoustic  boundary- 
layer  thickness.  In  a  latter  paper,  Culick3  connects  the  "flow-turning"  losses  to  an  increase 
in  entropy.  Since  the  publication  of  his  papers,  Flandro5,  Vuillot  and  Kuentzmann6  and 
Baum  and  Levine7  have  investigated  "flow-turning"  losses  within  the  acoustic  boundary- 
layers  associated  with  the  injection  of  fluid  into  a  chamber  containing  sound. 

Flandro,  using  an  admittance  correction  acoustic  "boundary-layer"  type  approach, 
extended  Culick's  longitudinal  model  to  include  the  acoustic  energy  losses  generated  by 
sound  waves  approaching  a  burning  wall  at  normal  as  well  as  grazing  flow  angles.  He 
derived  analytical  expressions  that  showed  that  the  acoustic  energy  losses  generated  from 
the  normally  directed  acoustic  waves  were  larger  than  the  the  losses  generated  from  the 
grazing  incident  sound  waves.  His  analysis  was  restricted  to  the  edge  of  the  viscous 
boundary-layer  within  which  the  mean  flow  direction  was  normal.  Thus  Flandro's  model 
assumed  all  the  acoustic  energy  losses  occurred  within  this  region  implying  that  other 
losses  associated  with  interactions  between,  say,  axial  mean  flow  gradients  and  the  sound 
field  within  the  interior  away  from  the  walls  were  apparently  negligible.  Flandro  estimated 
the  contribution  of  the  acoustic  losses  from  both  his  and  Culick's  model  to  T-Bumer 
stability  data  and  found  that  his  model  agreed  "better"  than  Culick's  model. 

Vuillot  and  Kuentzmann  conducted  a  series  of  experiments  to  examine  the  validity 
of  both  Culick's  and  Flandro's  models.  Test  results  favored  the  admittance  correction 
approach  adopted  by  Flandro.  In  a  separate  still  on-going  study,  Baum  and  Levine 
undertook  a  very  sophisticated  computational  investigation,  based  on  solutions  of  the  full 
Navier-Stokes  equations  of  motion,  of  the  effects  of  the  transverse  injection  of  air  into  a 
cylindrical  tube  containing  longitudinal  sound.  They  found  significant  sound  absorption 
takes  place  within  very  thin  regions  near  the  tube  wall  where  the  effects  of  viscosity  are 
important.  These  findings  support  Culick's  hypothesis  that  most  of  the  sound  absorption 
takes  place  within  very  thin  regions  near  the  wall.  They  also  found  (i)  the  sound 
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absorption  to  be  significantly  affected  by  refraction  induced  by  the  gradients  of  the  mean 
How,  (ii)  the  Richardson  effect  was  recovered  and  (iii)  acoustical  streaming  was  taking 
place  via  the  transfer  of  acoustic  energy  into  mean  flow  energy.  To  simplify  the  boundary 
conditions,  the  authors  assumed  the  tube  walls,  through  which  fluid  was  injected,  to  be 
acoustically  rigid.  Since  finite  admittance  side  walls  are  known  to  locally  affect  the  sound 
field,  the  effect  of  this  assumption  is  not  known. 

In  support  of  Culick's  rigid  wall  assumption,  it  is  only  proper  to  point-out  that  most 
of  the  researchers  in  this  field,  with  the  exception  of  Flandro,  assumed  that  the  injecting 
surfaces  were  acoustically  rigid.  In  a  real  injection  surface,  however,  the  local  wall 
admittances  are  finite  and  certainly  affect  the  behavior  of  the  sound  fields.  For  example,  if 
fluid  is  injected  through  small  orifices,  then  the  admittance  of  each  orifice  completely 
controls  the  behavior  of  the  local  sound  field  within  a  small  region  on  the  order  of  2-3 
orifice  diameters.  Unless  the  orifice  diameters  are  extremely  small,  these  regions  extend 
beyond  the  local  acoustic  boundary-layer  thickness  and  clearly  must  be  accounted  for  in  a 
theoretical  model. 

Culick2  applied  his  one-dimensional  "flow-turning"  loss  ideas  to  three-dimensional 
models.  He  claims  that  "owing  to  the  turning  of  the  flow  near  the  boundary,  there  are 
contributions  to  the  stability  of  waves,  which  do  not  arise  in  the  inviscid  three-dimensional 
analyses,  but  which  can  be  estimated  with  the  one-dimensional  analysis".  Culick 
incorporated  into  a  three-dimensional  analysis  parts  of  his  earlier  derived  one-dimensional 
"flow-turning"  loss  analysis.  This  was  accomplished  by  adding  additional  one¬ 
dimensional  terms  to  the  three-dimensional  formulas  for  the  complex  wave-numbers 
describing  the  two  classes  of  problems.  From  his  analysis  of  the  equations  describing  one¬ 
dimensional  complex  wave-numbers,  Culick  derived  analytic  expressions  for  the  acoustic 
velocity  near  a  boundary.  He  modelled  the  general  behavior  of  an  acoustic  velocity  incident 
to  a  surface  by  partitioning  it  into  parallel  and  non  parallel  components. 

Apparently  motivated  by  the  above  approach  taken  by  Culick,  Van  Moorhem8 
developed  a  procedure  to  allow  comparison  of  Culick’s  one  dimensional  wave  equation  to 
one  he  derived  by  averaging  the  three-dimensional  wave  equation  over  the  cross-section  of 
a  motor  cavity.  Applying  this  procedure.  Van  Moorhem  found  the  one-dimensional  model 
dws  not  contain  any  physical  phenomena  missing  from  the  three-dimensional  model  and 
that  the  addition  of  one-dimensional  "flow-turning"  terms  to  the  three-dimensional  wave 
equation,  as  claimed  by  Culick,  is  inappropriate.  Based  on  his  findings,  Van  Moorhem 
claims  that  Culick's  efforts  to  incorporate  some  of  his  "flow-turning"  terms  into  a  three- 
dimensional  stability  model  are  not  warranted 

It  is  not  clear  why  most  of  the  above  researchers  restricted  their  modelling  and 
prediction  of  "flow-turning"  losses  to  the  acoustic  boundary-layer  region.  Over  twenty 
years  ago,  Culick9  showed  that  the  streamlines  generated  from  the  injection  of  fluid  into  a 
rocket  motor  penetrates  into  its  center.  If  sound  absorption  takes  place  as  a  consequence  of 
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the  "flow-turning"  process,  then  it  appears  reasonable  to  assume  that  a  significant  fraction 
of  the  energy  loss  takes  place  within  the  rocket  interior  and  not  within  the  acoustic 
boundary  layer. 

In  an  effort  to  substantiate  Culick's  model,  Hersh  and  Walker10  undertook  an 
experimental  program  to  measure  "flow-turning"  losses.  The  results  of  their  study  showed 
that  sound  was  indeed  absorbed  by  the  transverse  injection  of  air  into  a  duct  containing 
longitudinal  sound.  Test  data  showed  that  part  of  the  sound  absorption  takes  place  at  the 
finite  admittance  side  walls  through  which  the  air  was  injected  and  part  within  the  duct 
interior.  Hersh  and  Walker  identified  refraction  of  the  sound  pressure  distribution  across 
the  duct  test  section  by  the  injected  flow  velocity  gradients  as  a  possible  important 
contributor  to  the  overall  sound  absorption.  The  idea  here  was  that  the  refraction  process 
increased  the  amplitude  of  the  sound  pressure  near  the  wall  thereby  increasing  the  local 
sound  absorption  by  the  side-walls. 

Hersh  and  Walker's  investigation  encountered  experimental  difficulties  that  may 
have  compromised  some  of  their  data.  In  order  to  generate  high  injection  flow  rates  into 
the  test  section,  two  compressors  were  connected  in  series.  This  caused  serious 
compressive  hearing  problems  requiring  the  installation  of  cooling  jackets  to  prevent  the  air 
temperatures  from  exceeding  the  safe  operating  limits  of  the  B&K  acoustic  microphones. 
At  the  higher  injection  flow  rates,  however,  the  fluid  temperatures  were  still  sufficiently 
close  to  the  maximum  allowable  values  to  cause  concern. 

Although  almost  the  same  experimental  facilities  used  by  Hersh  and  Walker  are 
used  in  this  investigation,  the  above  experimental  difficulties  were  avoided  by  using  suction 
rather  than  blowing  to  force  air  into  the  test  section.  This  eliminated  the  severe 
compressive  heating  of  the  air  as  it  was  injected  into  the  test  section.  The  experimental 
facilities  and  data  acquisition  and  reduction  scheme  used  in  this  investigation  are  described 
in  Section  2.  The  "flow-turning”  test  results  are  presented  in  Section  3.  The  test  results  are 
compared  to  a  one-dimensional  theoretical  model  derived  in  Section  4.  The  important 
findings  of  this  investigation  are  summarized  in  Section  5.  The  report  closes  with  a 
detailed  derivation  of  Culick’s  "flow-turning"  loss  expression  in  the  Appendix. 

EXPERIMENTAL  PROGRAM 


The  objective  of  the  experimental  program  was  to  measure  the  "flow-turning" 
losses  as  a  function  of  injection  flow  rate  and  sound  frequency.  According  to  Culick’s 
model,  the  sound  absorption  should  (i)  increase  almost  linearly  with  injection  flow  rate  at 
low  injection  rates  and  (ii)  be  independent  of  sound  frequency  because  the  interaction  takes 
place  between  the  injected  flow  and  longitudinal  (plane-wave)  sound.  The  experimental 
program  consisted  of  measuring  the  acoustic  energy  flux  across  a  test  section  as  a  function 


of  injection  flow  rate  for  test  frequencies  of  500  Hz,  1  KHz  and  2  KHz.  If  Culick’s  model 
is  valid,  than  sound  energy  should  be  absorbed  within  the  injection  test  section.  The 
acoustic  energy  flux  downstream  of  the  injection  test  section  should  be  less  than  its 
corresponding  upstream  value.  Further,  the  amount  of  energy  reduction  should  be  (1) 
proportional  to  the  injection  flow  rate  and  (2)  independent  of  the  sound  frequency. 

The  HAE  two-dimensional  "flow-turning''  facility  was  used  to  conduct  the 
measurements.  The  facility,  instrumentation  and  data  acquisition  scheme  are  described 
below. 

Facility  Description 

Figure  1  shows  a  scale  diagram  of  the  ''flow-turning"  facility.  The  test  section 
consists  of  a  3-inch  high  by  2.5-inch  wide  rectangular  cross-section  duct  with  sidewalls 
constructed  of  high  density,  1.25-inch  thick  Owens-Coming  Glasstrate  Fiberglass  board. 
This  material  approximates  reasonably  well  the  impedance  of  solid  motor  propellants.  Air 
flow  was  supplied  from  two  Paxton  CB-95  compressors  connected  in  series.  The 
compressors  were  operated  in  a  suction  mode  to  avoid  the  compressive  heating  problems 
encountered  in  the  earlier  study  by  Hersh  and  Walker.  The  air  flow  was  initially  sucked 
into  a  plenum  chamber  located  immediately  below  the  test  section.  Exiting  from  the 
plenum,  the  airflow  entered  the  test  section  via  the  eight  outlets  (four  on  each  side)  shown 
in  Figure  1.  Each  side  of  the  test  section  was  divided  into  four  segments,  with  individual 
air  plena,  control  valves  and  flow  monitoring  meters  allowing  uniform  transverse  flow  to 
be  injected  over  the  entire  wetted  area  of  the  Glasstrate  material.  This  was  necessary  to 
compare  test  results  to  a  theoretical  model  which  assumed  uniform  transverse  injection  of 
fluid  into  the  test  section. 

Upstream  of  the  porous- walled  test  section,  a  3-foot  long  hard  duct  was  connected 
to  the  JBL  2481  loudspeaker  shown  in  Figure  1.  The  loudspeaker,  used  to  introduce 
sound  into  the  test  section,  was  driven  at  frequencies  of  500,  1000  and  2000  Hz  by 
coupling  an  audio  generator  to  a  McIntosh  100  Watt  Power  Supply.  Downstream  of  the 
test  section,  a  6-foot  long  duct  was  connected  to  the  exponential  expansion  and  dissipative 
muffler.  The  muffler  was  used  to  reduce  the  amplitude  of  the  reflected  waves  as  much  as 
possible. 

Both  acoustic  pressure  and  hot-wire  velocity  measurements  were  recorded.  The 
acoustic  pressures  were  measured  using  a  B&K  4170  probe  tube  microphone  for  detailed 
surveys  of  the  sound  pressure.  This  instrument  combines  a  horn  coupled  microphone  with 
a  self-contained  preamplifier  for  measurements  in  confined  spaces.  This  instrument  was 
calibrated  against  a  B&K  4134  one-half  inch  microphone.  Entry  holes  for  the  probe 
microphone  were  drilled  at  appropriate  locations  upstream  and  downstream  of  the  test 
section,  allowing  three-dimensional  sampling  of  the  standing  wave  patterns  established 
within  these  sections. 
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The  acoustic  particle  velocities  were  measured  using  a  TSI  1280  hot-wire  probe 
connected  to  a  TSI  1054B  linearized  constant  temperature  anemometer.  The  hot-wire  was 
calibrated  against  a  pitot  tube. 

Data  Acquisition  Scheme 

The  data  was  collected  and  processed  using  a  MassComp  5400  computer  utilizing 
16  A/D  channels  and  12  D/A  channels  having  12-bit  resolution  integrated  into  the  computer 
to  perform  real  time  signal  processing  of  data.  MassComp  developed  software,  Laboratory 
WorkBench,  was  used  to  record  and  analyze  the  data.  A  block  diagram  of  the  data 
aquisition  scheme  is  shown  in  Figure  2.  The  A/D  converter  receives  two  signals,  one  of 
which  consists  of  the  output  from  a  sensor  [either  the  B&K  (pressure)  or  TSI  (velocity) 
instrumentation]  and  the  other  from  the  oscillator  driving  the  JBL  speaker  (this  output  is  not 
displayed  in  Figure  2  ).  The  signal  from  the  sensor  is  phase-averaged  to  remove  random 
noise  generated  from  the  turbulent  flow  past  the  sensor.  A  FFT  of  the  resulting  phase- 
averaged  signal  is  then  obtained  and  the  resulting  harmonic  pressure  component  recorded  in 
files  for  later  processing  and  analysis. 


TEST  RESULTS 


The  test  program  consisted  of  measuring  the  acoustic  energy  flux  across  the 
injection  test  section  as  a  function  of  injection  speed  and  sound  frequency.  Only  plane- 
wave  sound  was  generated  to  simplify  interpretation  of  the  data.  The  test  program  was 
divided  into  four  parts.  The  first  part  consisted  of  the  measurement  of  the  mean  velocity 
profiles  downstream  of  the  injection  test  section  as  a  function  of  injection  flow  speed. 
These  profiles  were  used  as  input  to  the  theoretical  model  presented  in  Section  4.  The 
second  part  consisted  of  the  measurement  of  the  acoustic  pressure  profiles  upstream  and 
downstream  of  the  injection  test  section  as  a  function  of  injection  speed  and  sound 
frequency.  These  measurements  were  used  in  the  energy  flux  calculations.  They  were  also 
used  to  determine  the  extent  the  sound  pressure  was  redistributed  across  the  test  section  by 
refraction  by  the  injected  flow  velocity  gradients.  The  third  part  consisted  of  the 
measurement  of  the  acoustic  velocity  profiles  downstream  of  the  injection  test  section  as  a 
function  of  injection  flow  speed  and  sound  frequency.  These  measurements  were  used  as 
an  independent  check  of  the  pressure  refraction  findings.  The  fourth  and  final  part 
consisted  of  measuring  the  standing  wave  pattern  upstream  and  downstream  of  the  injection 
test  section.  These  measurements  were  used  to  calculate  acoustic  energy  flux.  The  results 
of  these  tests  are  presented  below. 
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Linearized 


Mean  Velocity  Profiles 

The  first  task,  undertaken  was  to  survey  the  mean  velocity  profiles  downstream  of 
the  test  section  as  a  function  of  injection  Mach  number.  A  TSI-1054B  linearized  constant 
temperature  anemometer  was  used  to  measure  the  velocity  profiles.  It  was  calibrated 
against  a  standard  pitot  tube  before  each  test  run  to  ensure  its  linear  response  against  flow 
speed.  A  typical  calibration  curve,  displayed  in  Figure  3,  shows  a  linear  calibration  curve 
for  flow  speeds  beyond  about  3  m/s.  Below  this  limit,  however,  the  slope  of  the 
calibration  curve  is  slightly  reduced.  This  is  probably  caused  by  thermal  conduction  effects 
at  low  speeds  within  the  facility. 


Figure  3.  Calibration  of  Hot-Wire 


The  mean  velocity  profiles  immediately  downstream  of  the  injection  test  section  are 
displayed  in  Figure  4  as  a  function  of  injection  flow  speed.  Observe  that  the  flow  is  highly 
rotational  with  very  thick  regions  of  shear.  The  magnitude  of  the  shear  increases  with 
increasing  injection  flow  rate. 

Acoustic  Pressure  Profiles 

Measurements  of  the  sound  pressure  distributions  across  the  duct  centerline  - 
upstream  and  downstream  of  the  injection  test  section  -  were  undertaken  to  verify  that  only 
plane-wave  sound  was  generated  at  the  test  frequencies  of  500  Hz,  I  KHz  and  2  KHz. 
The  B&K  Type  4170  Probe  Microphone  described  in  Section  2  was  used  to  measure  the 
sound  pressure  amplitude. 
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Figure  4.  Downstream  Mean  Velocity  Profiles 


The  results  of  the  survey,  upstream  of  the  test  section,  are  displayed  in  Figure  5. 
These  measurements  were  conducted  in  the  absence  of  injected  flow.  The  acoustic 
pressure  distributions  across  the  duct  at  three  different  heights  were  recorded.  The  data 
show  that  plane-wave  sound  was  generated.  These  results  were  anticipated  because  the 
frequencies  were  below  the  duct  first  transverse  mode  cut-on  eigenfrequency. 

The  corresponding  sound  pressure  distributions  downstream  of  the  test  section 
were  measured  for  injection  Mach  number  of  0.001,  0.002  and  0.003.  The  interest  here 
was  whether  the  mean  injected  flow  velocity  gradients  redistribute  the  sound  pressure 
amplitudes  across  the  test  section.  The  redistribution  of  sound  pressure  amplitude  can  be 
explained  in  terms  of  refraction  by  the  mean  flow  velocity  gradients.  Figure  6  illustrates 
the  velocity  refraction  mechanism.  For  the  case  of  downstream  sound  propagation,  the 
phase  speed  is  greater  in  the  center  than  near  the  walls  due  to  the  mean  flow  profile.  The 
higher  phase  speed  at  the  center  turns  the  wave  front  into  the  wall  thereby  decreasing  the 
sound  pressure  at  the  duct  center  and  increasing  it  at  the  wall.  The  converse  is  true  for  the 
case  of  upstream  sound  propagation.  Here  the  sound  pressure  is  greater  in  the  duct  center 
than  at  the  wall. 
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Hersh  and  Walker  showed  that  the  response  of  a  probe  tube  microphone  was 
strongly  affected  by  the  grazing  flow  past  its  tip.  The  grazing  flow  increases  the  acoustic 
resistance  of  the  microphone  in  the  same  way  that  the  resistance  of  a  Helmholtz  resonator  is 
increased  by  a  grazing  flow.  This  increase  in  resistance  results  in  a  decrease  in  the 
response  of  the  probe  tube.  The  influence  of  the  grazing  flow  was  corrected  by  calibrating 
the  B&K-4170  probe  tube  against  a  B&K-4134  one-half  inch  microphone.  This  was 
accomplished  by  flush  mounting  the  4134  microphone  in  the  side  wall  of  the  HAE  Wind 
Tunnel  and  positioning  the  4170  probe  tube  in  the  air  stream  immediately  in  front  of  the 
4134  as  shown  schematically  in  Figure  7.  The  response  of  the  microphones  was  measured 
as  a  function  of  injection  flow  speed  at  test  frequencies  of  500  Hz,  1  KHz  and  2  KHz.  The 
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results  of  this  calibration  procedure  are  displayed  in  Figures  8  &  9.  Figure  8  displays  the 
ratio  of  flush  mounted-to-probe  tube  sound  pressure  amplitude  as  a  function  of  grazing 
flow  speed  for  the  three  test  frequencies.  The  data  have  been  normalized  by  the  zero  speed 
values  of  each  microphone  and  show  that  only  modest  corrections  are  necessary,  the  largest 
being  about  9%  for  the  1  KHz  data.  Figure  9  shows  the  corresponding  phase  corrections 
representing  the  differences  in  phase  between  the  4134  and  4170  microphones.  As  shown, 
these  corrections  are  also  modest. 


B&K  4170  Probe  Tube 
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Preliminary  measurements  showed  the  4170  probe  tube  to  be  sensitive  to  mean 
flow  speed  and  mean  flow  shear.  To  avoid  extensive  calibration  of  the  probe  tube,  the 
B&K  UA-0385  bullet  shaped  nose  cone  was  attached  to  a  B&K-4136  one-quarter  inch 
microphone  to  measure  the  sound  pressure  profiles  downstream  of  the  test  section.  As 
displayed  in  Figures  10  &  11,  only  plane-wave  sound  was  generated  across  the  test 
section,  even  for  the  maximum  injection  Mach  number  tested  of  0.003.  Thus  the  effects  of 
refraction  by  the  mean  flow  velocity  gradients  were  negligible  -  at  least  for  the  range  of 
injection  flow  speeds  tested. 

The  acoustic  pressure  measurements  presented  above  suggest  that  plane-wave 
motion  exists  throughout  the  test  section.  Additional  sound  pressure  surveys  were 
conducted  in  the  vertical  direction  (see  Figure  1)  to  verify  that  the  sound  field  is  planar  over 
the  entire  duct  cross-section.  The  B&K  4170  probe  tube  was  used  to  measure  the 
amplitude  and  phase  distributions  in  the  vertical  direction  within  the  injection  test  section 
exit  plane.  The  results  of  the  pressure  sun-ey,  presented  in  Figures  12  &  13  for  the  three 
test  frequencies,  show  that  the  acoustic  motion  throughout  the  entire  injection  test  section  is 
planar. 

Acoustic  Velocity  Profiles 

The  TSI  hot-wire  anemometry  system  was  used  to  measure  the  acoustic  particle 
velocity  profile  across  the  test  section  to  provide  an  independent  verification  that  the  sound 
wave  refraction  by  the  mean  flow  velocity  gradients  was  indeed  negligible.  The  idea  here 
is  that  for  sound  fields  governed  by  linear  acoustics,  a  linear  relationship  exists  between  the 
sound  pressure  p'  and  the  particle  velocity  u\  i.e.  p'  ~  u’.  Thus,  both  the  acoustic 
pressure  and  velocity  profiles  should  be  planar.  Preliminary  measurements  were  conducted 
to  verify  that  p'  ~  u'.  The  amplitudes  of  p'  and  u'  were  measured  as  a  function  of  speaker 
input  voltage,  the  results  of  which  are  plotted  in  Figure  14  for  frequencies  of  500  Hz,  1 
KHz  and  2  KHz.  Observe  that  both  p'  and  u'  increase  linearly  with  sound  pressure 
amplitude  which  is  proportional  to  speaker  input  voltage. 

The  linear  response  of  both  these  measurements  verified  that  the  hot-wire  measured 
acoustic  velocities.  The  acoustic  particle  velocity  profiles  across  the  test  section,  displayed 
in  Figures  15  &  16,  show  that  the  sound  is  plane-wave  in  support  of  the  findings  displayed 
in  Figures  10  &  11.  The  data  near  the  wall  is  not  accurate  because  the  hot-wire  output  is 
not  accurate  when  the  mean  flow  velocity  is  small. 

Acoustic  Energy  Transmission  Measurements 

According  to  Culick's  "flow-turning"  model,  the  interaction  between  transversely 
injected  mean  flow  and  longitudinally  excited  plane-wave  sound  results  in  a  net  absorption 
of  sound.  To  substantiate  his  findings,  an  acoustic  energy  exchange  concept  (AEE)  is 
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Figure  10,  Downstream  Sound  Pressure  Amplitude  Survey 


Figure  11.  Downstream  Sound  Pressure  Phase  Survey 
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introduced  below  defined  as  the  difference  between  the  acoustic  energy  transmission 
coefficient  (T)  across  the  test  section  with  and  without  injected  mean  flow, 


AEE  =  T(Mmj)  -T(0)  (2) 

Here  T(Mjnj)  is  the  acoustic  energy  transmission  coefficient  across  the  test  section  defined 
as 


T(M|ni )  ,  lOLosj 


P.(C  M,„,  )u„(L,  M,„) 
P»(0.  Mi,,)ur(°.  m,„,) 


(1  +  M„) 


(3) 


where  pr  and  ur  represent  the  real  parts  of  the  root-mean-square  (mas)  acoustic  pressure 
and  velocity  respectively  and  Mav  is  the  average  Mach  number  across  the  duct  defined  as 
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The  parameter  H  is  the  duct  width  and  M(y)  is  the  mean  flow  Mach  number  distribution 
across  the  duct  cross-section  displayed  in  Figure  3.  The  term  1+Mav  in  Eq.(3)  accounts  for 
the  additional  energy  convection  of  the  sound  field  by  the  mean  flow11. 

Since  the  data  presented  in  Section  2.3  showed  that  the  sound  field  across  the 
injection  test  section  was  planar,  the  acoustic  pressures  and  velocities  are  evaluated  at  the 
duct  center  (y=32  mm).  This  simplifies  enormously  the  time  and  effort  required  to  conduct 
the  acoustic  energy  transmission  loss  measurements.  The  rms  quantities  in  Eq.  (3)  are 
defined  in  the  usual  way  as 
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The  B&K  4170  probe  tube  was  used  to  measure  the  sound  pressure  standing  wave 
patterns  upstream  and  downstream  of  the  injection  test  section  for  sound  frequencies  of  500 
Hz,  1  KHz  and  2  KHz  and  for  injection  Mach  numbers  of  0.001,  0.002  and  0.003.  Test 
results,  uncorrected  for  the  effects  of  grazing  flow  past  the  probe  tube,  are  displayed  in 
Figures  17-22.  As  shown,  curve-fits  of  the  data  are  superimposed.  The  curve-fits  were 
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Figure  18.  Curve-fit  of  Standing  Wave  Pattern  for  500  Hz 
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Eiqure  21,  Curve-fit  of  Standing  Wave  Pattern  for  2  KHz 
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derived  starting  from  the  following  standard  one-dimensional  solution  of  the  sound 
pressure  containing  incident  and  reflected  components  in  a  moving  fluid, 


p(x,  t )  =  A 


'(“*  +  rnr-  +  *} 

e  ^  ;  +  Re  ^  “  > 


(6) 


where  A  is  the  amplitude  of  the  incident  sound  pressure,  R  is  the  ratio  of  the  amplitudes  of 
the  incident  and  reflected  waves,  <)>  is  the  phase  of  the  reflected  wave  and  k  =  to/c.  The 
corresponding  one-dimensional  sound  particle  velocity  is 


u(x,t) 


cot  + 


kx 
1  +  M 

m 


(7) 


The  following  rms  sound  pressure  (Pnns)  is  derived  from  Eqs.  (5  &  6), 


P  rmi  00  ~  A 


f  1  +  R  +  2R  cos 


2kx 


(8) 


Observe  from  Eq.  (8)  that  minimums  (Pmjn)  occur  whenever  2kx+<j>  =  7t  +27in,  and 
maximums  (Pmax)  occur  whenever  2kx+<(>  =  27tn  where  n  is  an  integer.  Using  this,  it  is 
straight-forward  to  show  that 


P  +  P 

min  ">« 


P  ~  P 

R  _  _  r  min 


2kx 


+  P 


=  7t  - 


min 


min 


1-M*„ 


(9) 


The  final  estimate  of  the  standing  wave  pressures  upstream  and  downstream  of  the 
injection  test  section  (see  Fig.  1)  follows  upon  substituting  Eq.  (9)  into  Eq.  (8).  As  shown 
in  Figures  17-22,  the  agreement  between  predicted  and  measured  standing  wave  patterns  is 
excellent.  Incorporating  the  above  definitions  into  Eq.  (3)  yields  the  following  expression 
for  the  acoustic  energy  transmission  T  across  the  injection  test  section, 
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T(M,„i)  =  10 '°e 


A*0(L,  M,n| )  ( I  —  Rp)  (1  +  M„ ) 


(10) 


The  acoustic  energy  transmission  depends  only  upon  the  amplitude  A  of  the 
incident  wave,  the  coefficient  R  of  the  reflected  wave  and  the  mean  flow  average  Mach 
number.  Here,  Ad  &  Au  represent  the  amplitudes  of  the  sound  pressures  downstream 
(x=L)  and  upstream  (x=0)  of  the  injection  test  section  respectively.  Table  1  summarizes  the 
parameters  A  and  R  for  the  twenty-four  test  cases.  The  downstream  amplitudes  Ad  are 
corrected  for  the  effects  of  the  grazing  flow  past  the  tip  of  the  B&K  4170  Probe  Tube. 
These  corrections  were  provided  by  the  amplitude  calibration  curves  shown  in  Figure  8. 
Substituting  these  corrections  along  with  the  other  parameters  displayed  in  Table  1  into  Eq. 
(10)  results  in  the  acoustic  energy  transmission  values  summarized  in  Figure  23  and  Table 
2. 

The  results  show,  for  the  three  frequencies  tested,  a  slight  acoustic  energy  gain  for 
injection  Mach  numbers  between  0  <  Minj  <0.1.  Beyond  Minj  =  0.1,  the  data  shows 
virtually  little  change  for  frequencies  of  500  Hz  and  1  KHz  but  over  1  dB  of  loss  at  2  KHz. 
This  is  more  clearly  defined  in  Figure  24  which  was  derived  using  the  acoustic  energy 
exchange  definition  of  Eq.  (2).  What  may  be  happening  here  is  that  at  very  low  injection 
Mach  numbers,  the  absorption  of  sound  energy  may  be  too  little  to  off-set  the  increase  in 
energy  due  to  mean  flow  convection  1+Mav  as  indicated  in  Eq.  (10).  This  may  explain  the 
slight  acoustic  energy  gain.  At  the  higher  injection  flow  Mach  numbers,  however,  the 
flow-turning"  losses  may  increase  more  rapidly  than  the  convection  term  resulting  in  the 
small  losses  shown  in  Figures  23  &  24.  If  this  explanation  is  valid,  it  suggests  that  the 
"flow-turning"  losses  are  not  controlled  by  terms  proportional  to  the  first  power  of  Minj,  as 
Culick's  linear  analysis  predicts,  but  instead  exhibits  a  quadratic  or  higher  power 
dependence  upon  Mjnj.  Regardless  of  the  details  of  the  sound  absorption  process,  it  seems 
reasonable  to  conclude  from  the  measurements  presented  in  Figures  23  &  24  that  the  "flow¬ 
turning"  process  indeed  absorbs  sound.  This  absorption  is  certainly  very  small  relative  to 
the  absorption  by  the  finite  admittance  side  walls  (at  Mjnj  =  0). 


THEORETICAL  MODEL 


The  development  of  a  detailed  one -dimensional  theoretical  model  was  undertaken  to 
extend  Culick's  rigid  wall  derivation  to  include  the  effects  of  finite  admittance  (i.e., 
absorbing)  walls.  The  interest  here  is  to  identify  and  rank  order  the  various  interactions 
among  the  injected  flows,  longitudinal  sound  fields  and  finite  admittance  walls.  Another 
equally  important  objective  is  to  examine  the  value  and  utility  of  the  one-dimensional  model 
by  comparing  its  predicted  absorption  to  the  measured  data  presented  in  Section  3. 
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Table  1.  Summary  of  Uncorrected  Standing  Wave  Pattern  Test  Results 


frequency  Minj  Upstream  Downstream 

(Hz)  A  R  Aunc  R 


500 

0 

36,344 

0.315 

27,100 

0.429 

0.001 

22,830 

0.335 

18,350 

0.513 

0.002 

27,800 

0.334 

22,211 

0.525 

0.003 

20,975 

0.310 

15,830 

0.539 

1  KHz 

0 

129,300 

0.115 

92,404 

0.296 

0.001 

88,622 

0.134 

64,520 

0.305 

0.002 

87,124 

0.137 

60,120 

0.329 

0.003 

72,493 

0.113 

46,750 

0.338 

2  KHz 

0 

34,551 

0.075 

21,170 

0.106 

0.001 

24,010 

0.073 

14,840 

0.106 

0.002 

24,992 

0.061 

13,889 

0.113 

0.003 

17,917 

0.072 

9,031 

0.129 

Table  2.  Summary  of  Transmission  Measurements 


Frequency 

(Hz) 

Minj 

Downstream 

Auncorr  Aeon- 

T(Minj) 

(dB) 

Acoustic  Energy 
Exchange  -  (dB) 

500 

0 

36,344 

27.100 

-2.98 

0.00 

0.001 

22,830 

18,405 

-2.62 

0.36 

0.002 

27,800 

22,544 

-2.60 

0.38 

0.003 

20,975 

16,432 

-3.02 

-0.04 

1  KHz 

0 

129,300 

92,404 

-3.26 

0.00 

0.001 

88,622 

66,843 

-2.73 

0.53 

0.002 

87,124 

64,449 

-2.92 

0.54 

0.003 

72,493 

52,313 

-3.16 

0.10 

2  KHz 

0 

34,551 

21,170 

-4.28 

0.00 

0.001 

24,010 

14,899 

-4.11 

0.17 

0.002 

24,992 

14,167 

-4.86 

-0.58 
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AEE(Minj;f)  -  dB 


Figure  24.  Measured  Acoustic  Energy  Exchange 
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Several  important  assumptions  were  introduced  into  the  model  to  simplify  its 
derivation.  First,  the  channel  into  which  fluid  is  injected  is  uniform  in  cross-section  and 
infinite  in  length.  This  assumption  was  required  because  Culick's  model  assumed  that  the 
rocket  chambers  were  bounded  by  rigid  walls  at  both  ends.  Second,  longitudinal  sound  is 
assumed  to  be  established  within  the  channel  prior  to  the  injection  of  the  fluid.  Since  the 
channel  is  infinite  in  length,  only  propagating  sound  waves  are  established  within  the 
channel.  Third,  the  fluid  is  inviscid  and  non-heat  conducting.  Fourth,  the  temperature  of 
the  injected  fluid  is  equal  to  the  fluid  temperature  within  the  channel.  Fifth,  the  perfect  gas 
law  is  used  as  the  equation  of  state  of  the  injected  fluid.  And  sixth,  fluid  is  injected  at  very 
low  volume  flow  rates  into  the  channel  in  a  uniform  manner  along  its  upper  and  lower 
w  alls.  The  fluid  is  injected  in  a  direction  perpendicular  to  the  channel  axis. 

Derivation  of  Governing  Acoustic  Equations 

A  schematic  of  the  control  volume  used  in  this  derivation  is  shown  in  Figure  25. 
Mean  flow  is  injected  into  the  control  volume  at  y  =  0  and  y  =  H  as  shown.  Assuming  the 
fluid  is  inviscid,  the  conservation  of  mass  and  axial  momentum  may  be  written, 
respectively,  as 


Injected  Mass  Flow  Rate 
pv 


Injected  Mass  Flow  Rate 


5p 

at 


+ 


^-(pu)  =  gpv 


and 


au 

at 


ap 

ax 


=  -Hpuv 


(11) 


(12) 


The  derivation  of  Eq.  (11)  assumes  that  mass  is  pumped  transversely,  as  shown 
schematically  in  Figure  25,  into  the  control  volume  at  the  rate  pv  per  unit  area.  Adopting 
the  usual  practice  of  assuming  that  the  instantaneous  fluid  properties  can  be  decomposed 
into  mean  and  acoustic  components,  the  following  conservation  of  acoustic  mass  flux  is 
derived, 


at 


T5T  +  ui£r +  ■jsrP' =  ■r(vp'  +  pv') 


au- 
+  p  ax 


(13) 


Here  (')  denotes  acoustic  quantities  and  p  denotes  the  constant  mean  fluid  density.  The 
corresponding  acoustic  momentum  flux  equation  is 


at 


u 


ay  au 

ax  ax 


u' 


+  =  ~H[pUv'  +  V(i’u'+Up')l 


(14) 


The  local  mean  flow  speed  [LfeU(x)]  at  a  distance  x  from  the  leading  edge  of  the 
injection  region  is  calculated  from  the  conservation  of  volume  flow.  This  is  valid  because 
of  the  very  low  rates  at  which  the  fluids  are  injected  into  the  channel.  Referring  to  Figure 
25,  the  total  volume  flow  rate  injected  into  the  channel  from  the  upper  and  lower  surfaces  is 
2xWV  where  V  is  the  constant  injection  speed  and  W  is  the  channel  width.  The  local  mean 
flow  speed  U(x)  is  calculated  form  the  conservation  of  mean  volume  flow 


2xWV  =  WHUx)  ->  IKx)  =  2  £  V 

H 


(15) 


It  is  of  interest  to  observe  that  this  approach  agrees  exactly  with  the  more  accurate 
two-dimensional  analysis  developed  by  Culick9.  The  derivation  of  his  model  starts  with 
the  assumption  that  the  injected  mean  flow  is  two-dimensional,  inviscid  and 
incompressible.  The  governing  mean  flow  continuity  equation  is 
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(16) 


dU  ,  dV 
dx  dy 


0 


Adopting  the  approach  used  by  Culick ,  the  following  solution  was  found. 


U(x.Y)  =  Vtt  (**->in  (2L);  V  (y )  =  V„  cos(^-)  (17a,b) 


The  accuracy  of  Eq.  (17a)  was  verified  by  comparing  it  to  measured  values  of  the 
mean  flow  profile,  downstream  of  the  injection  test  section  where  x  =  16  inches.  As 
shown  in  Figure  26,  the  agreement  between  prediction  and  measurement  is  excellent 
everywhere  except  near  the  walls  where  the  assumption  that  the  mean  flow  is  inviscid  is  not 
expected  to  be  valid.  Observe  that  the  predicted  axial  velocity  U(x,y)  satisfies  the  no-slip 
condition  at  the  walls  and  that  V(y)  satisfies  the  inward  injection  speeds  at  y  =  0  and  y  =  H. 
Observe  further  that  the  successful  prediction  of  the  outer  portion  of  the  velocity  profiles 
suggests  that  the  injected  mean  flow  streamline  patterns  are  laminar.  The  average  speed 
across  the  channel  is  calculated  by  substituting  Eq.(17a)  into  Eq.  (4)  to  yield 


Figure  26.  Plots  of  Predicted  &  Measured  Velocity  Profiles 
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(18) 


* 


which  agrees  exactly  with  Eq.  (15). 

Since  the  fluid  temperatures  inside  and  outside  of  the  test  facility  are  assumed  to  be 
equal,  the  conservation  of  energy  flux  may  be  written 


+ 


£(PU)  =  fjPv(h 


(19) 


Here  e  =  CVT  is  the  fluid  internal  energy  per  unit  mass  and  h  =  CpT  is  the  fluid  internal 
enthalpy  per  unit  mass.  Replacing  T  by  the  equation  of  state  for  ideal  gases. 


P  =  pRT  = 


C2  =yRT 


(20a, b) 


the  following  conservation  of  energy  flux  may  be  derived  after  some  algebra, 


Expanding  Eq.  (21)  into  mean  and  acoustic  components  and  retaining  terms  linear 
in  mean  flow  speed  U  yields  the  following  acoustic  conservation  of  energy  flux. 


Perturbation  Solution 

Observe  that  Eqs.  (13,  14  &  22)  contain  the  four  unknowns,  p',p',u'  and  v'.  A 
fourth  equation  must  be  introduced  in  order  to  solve  these  equations.  Since  a  one- 
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dimensional  solution  is  sought,  v'  is  connected  to  the  sound  pressure  p'  through  the  non- 
dimensionalized  admittance  of  the  side  wall  panels  Aw  =  Ar  +  iAj,  as  defined  below, 


v'  = 


-(A„+iA,)£ 


(23) 


Referring  to  Figure  25,  in  order  for  sound  energy  to  be  absorbed,  the  local  acoustic 
velocity  must  be  pumped  into  the  walls  at  y  =  0  and  y  =  H.  Thus  v’  must  be  negative  as 
shown  in  Eq.  (23).  Replacing  v'  by  Eq.  (23)  and  substituting  into  Eqs.  (13  &  14)  yields 
the  following  equations  governing  the  conservation  of  acoustic  mass  and  axial  momentum, 


at 


and 


ap* 

ax 


(24) 


(25) 


Equations  (24  &  25)  can  be  combined  into  a  single  higher-order  equation  by  differentiating 
Eq.  (24)  by  d/dt  and  Eq.  (25)  by  d/dx  yielding 


aV 

at2 


aV 

ax2 


2A  w  ap' 
ch  at 


2V 

Aw 

a2u*  ay 

aw  ap'" 

CH 

[pc3ir- 

"2_hp'J  +  > 

i  pc — -  -  c — — 

.  ax2  axat 

1  h  ax 

(26) 


Equation  (24)  is  used  to  eliminate  u'  from  Eq.  (26)  resulting  in  two  equations  governing 
the  time  and  spatial  changes  of  p'  and  p'. 

Following  the  approaches  adopted  by  Culick  and  Flandro,  approximate  solutions  to 
Eqs.  (22  &  26)  are  sought  based  on  the  following  regular  power  series  expansion  of  p'  and 
p’  in  terms  of  the  small  parameter  Mi„j, 


t 
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and 


p'(x,  t)  2  p,(x,  t)  +  Mln,  p, (x,  t)  +  Mjn, p2(x,  t)  +  ... 


(27) 


p'(x,  t)  =  pQ(x,  t)+  M  p^(x,  t)+  M  p2(x,  t)  + 


where 


Wx)  =  ^  =  M„(^);  M,„,  =  £ 


(28) 


(29a, b) 


Observe  that  the  (  )'  notation  used  to  denote  acoustic  variables  has  been  omitted  for 
convenience. 

Combining  the  resulting  equations  to  eliminate  p'  and  collecting  the  coefficients  of 
the  first  two  powers  of  M,nj  leads  to  a  collection  of  conservation  equations  modelling  the 
behavior  the  sound  field.  For  the  lower-order  terms  (independent  of  Mjnj),  these  equations 
are 


a 


=  o 


d'p.  _  ^P.  2A*  ap,  _ 

dx2  ax2  ch  at 


(30) 

(31) 


For  the  terms  proportional  to  the  corresponding  conservation  equations  are 


ajPi 

at 


(32) 


and 
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8aPl  8*pt  2Awapt 
at1  3x2  cH 


The  term  labelled  F-T  in  Eq.  (33)  represents  the  "flow-turning"  term  used  in  the 
Appendix  to  derive  Eq.  (1).  This  term  will  later  be  shown  to  dominate  the  behavior  of  the 
interaction  between  the  injected  fluid  and  the  longitudinal  sound  field  within  the  channel. 

Equations  (30  &  32)  show  that  the  acoustic  pressure  and  density  of  the  lowest- 
ordered  and  linear  terms  are  adiabatically  related  where 

p#  =  e2p#  and  pt=c2pt  (34a,b) 


Substituting  Eq.  (34a)  into  Eq.  (31)  yields  the  following  lower-order  wave  equation 
governing  the  acoustic  behavior  of  the  sound  pressure 


1  d*P.  dP.  2AW  3P. 

c2  at2  3x2  31 


(35) 


The  corresponding  first-order  wave  equation  is  determined  by  substituting  Eq.  (34b)  into 
Eq.  (33)  to  yield 


1  d*Pi  a2P,  2AW  gPi 
C2  at2  3x2  CH  at 


(36) 
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An  approximate  solution  to  Eq.  (35),  valid  for  small  values  of  Ar  and  Ap  is 


-  A 

P,(x,  t)  =C,e 


£ 

RH 


l(<o«  -  kx ) 


(37) 


w  here  Co  is  an  arbitrary  constant.  Equation  (37)  represents  an  approximate  solution  of  the 
sound  pressure  distribution  within  the  injection  test  section  in  the  absence  of  injection  flow. 
Referring  to  the  discussion  following  Eq.  (23)  ,  the  negative  sign  was  inserted  in  order  to 
force  the  solution  defined  by  Eq.  (37)  to  agree  with  the  measured  acoustic  energy 
transmission  data  (  at  M;nj  =  0)  shown  in  Figure  23.  No  attempt  has  been  made  in  the 
above  derivation  to  account  for  the  effects  of  reflected  waves.  Although  such  a  solution 
can  be  found,  the  algebra  is  very  messy  and  little  will  be  served.  Referring  to  Table  1,  the 
upstream  propagating  acoustic  energy,  which  is  proportional  to  the  square  of  the  reflection 
coefficient  R,  is  quite  small  for  ail  three  test  frequencies. 

The  acoustic  energy  transmission  ratios  across  the  test  section  in  the  absence  of 
injected  flow  can  be  easily  estimated  from  Eq.  (37).  Following  standard  acoustic  analysis, 
the  real  part  of  the  rms  sound  pressure  is 


p«(x,=  VTe 


»  L 


(38) 


The  acoustic  energy  transmission  coefficient  (T)  across  the  test  section  follows  from  Eqs. 
( 1 0  &  38)  to  be 


T(0)  =  10  Lo 


r  _2 

Pr 

9  -7 


(L,  0) 


LPr(0’0)J 


=  -8.  69^ A  R 


(39) 


Equation  (39)  is  used  in  conjunction  with  the  measured  T(0)  data  summarized  in  Table  2  to 
back-out  predicted  values  of  Ar.  These  values  are  compared  to  the  values  measured  by 
Hersh  and  Walker  in  Figure  27.  The  agreement  is  reasonable.  Since  measured  data  was 
used,  the  values  of  Ar  shown  in  Figure  27  include  the  effects  of  viscous  damping. 
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Figure  27.  Predicted  vs  Measured  Re(Admittance) 


The  solution  to  pt  is  obtained  from  Eq.  (36)  by  substituting  Eq.  (37)  forp# ,  Eq. 
(34a, b)  to  connect  p'  and  p’,  and  using  the  lowest-order  mass  flux  conservation  equation 
to  derive  the  following  connection  between  the  lowest-order  particle  velocity  uG  to  the 
pressure  pQ, 


1  dpa  du0 
c*  dx  P  dx 


2A„ 


p^pw^-dr(2A-  +  'kH  )p° 


cH 


(40) 


This  yields 


1  *  Pi  *  Pi  2A«  * i 
— —  + 


2  2 
C  dt 


dx 


cH  dt 


2Pn 


8  Ap  +  i(2kH  +  8/4,) 


■4Tf(kH  +  AMkH+2A>)’2iA4  (40 
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A  solution  to  Eq.  (41)  is  sought  of  the  form 


p;(x,  y,  t)  = 


+  i(<0»  -lot) 

f  (x) 


(42) 


Substituting  Eqs.  (37  &  42)  into  Eq.  (41)  and  assuming  that  both  Ar&Ai  «  1  leads  to  the 
following  simplified  differential  equation  governing  f(x), 


^{16A„  +  6ikH  +  4(kx)  [kH  -  2iA„  ]} 


Since  AR  «  1 ,  an  approximate  solution  to  Eq.  (43)  is  sought  in  terms  of  the  expansion 


f  (x)~f0(x)+A^(x)  +  A*f2(x)  +  „.. 


(44) 


Substitution  of  Eq.  (44)  into  Eq.  (43)  yields  the  following  differential  equations  to  order 
(Ar)2, 


f",  -2ikf'#=-  ~ /6ikH  +  4(kx)kHl 

H  ITT  / 


(45) 


f",  -  2ikf  'j  =  ~f\  +  -i(i6  -  8ikx) 

H 


(46) 


and 


-“'•■ft,.- 


H 


1*. 


(47) 
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The  "flow-turning"  term  appears  as  a  source  term  in  Eq.  (45).  Foregoing  the  details,  the 
following  particular  solution  was  derived 


The  composite  solution,  to  order  Mjnj,  of  the  one-dimensional  sound  pressure 
distribution  within  the  test  section  follows  by  substituting  Eqs  (37,  42  &  48)  into  Eq.  (27) 
to  yield 


p(x,t)  =C,e 


A  —  +  l(<o»  -  kx) 

"H 


[l+MInlf(x)] 


(49) 


The  root-mean-square  pressure  distribution  within  the  injection  test  section  follows 
from  application  of  Eq.  (5)  to  be 


P«(x’  M»,  )  =  “Y6  *"L[I+  M„,f  R(x)]  (50) 

where 


f  R(x)  -  2  Ar 

k  ■  ..  -y  i 

F-  T 

Here  fR(x)  represents  the  real  part  of  f(x).  The  solution  defined  by  Eq.  (51)  shows  that  the 
"flow-turning"  contribution  dominates  the  other  terms  providing  AR  «  1.  Further,  the 
negative  sign  indicates  that  sound  energy  is  absorbed.  Since  we  are  dealing  only  with 
travelling  waves,  the  acoustic  energy  transmission  across  the  test  section  is  expressed  as 


2(kH)  J 


+  A 


13 


4(kH) 


(51) 
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TfM , )  =  lOLog  -P - ^-(1  +  M3 


(52) 


Substituting  Eq.  (51)  into  Eq.  (50)  yields  the  following  expression  for  the  acoustic 
energy  transmission  across  the  test  section. 


T  (Mmj )  s  T  (0)  +  101.og  fl+  +10Log 


p  +  *V»(t>l 


(53) 


Rewriting  Eq.  (53)  using  Eq.  (2)  results  in  the  following  theoretical  estimate  of  the  acoustic 
energy  exchange  (AEE)  between  transversely  injected  flow  and  longitudinal  sound. 


AEE=  lOLog(l  +  M„v)  +  lOLog 


Equation  (54)  is  the  desired  result.  Its  interpretation  requires  considerable 
comment.  To  begin  with,  observe  that  since  AR  «  1,  Eq.  (54)  is  dominated  by  the  "flow- 
turning"  term  -2(L/H)Minj,  This  term  can  be  traced  to  Culick's  "flow-turning"  loss 
expression  defined  by  Eq.  (1)  for  the  special  case  of  rigid  walls  where  Ar  =0.  Under  this 
condition,  Eq.  (54)  simplifies  to 


AEE=  I0Log(l  +  Mav)  +  lOLo 


(55) 


Combining  Eqs.  (15  &  29a),  the  "flow-turning"  term  2Mjnj(L/H)  in  Equation  (55)  can  be 
simplified  by  observing  that 
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(56) 


Now  substituting  Eq.  (56)  into  Eq.(55)  produces  the  interesting  result  that 


AEE=  10Log/l  +  M„  +  10Log,l  -  M„  )  =  IOLog(  I  -  M'„v ) 

"  tttJ  <57) 


Convection 


Equation  (57)  suggests  that  within  the  context  of  Culick's  "flow-turning"  model, 
the  interaction  between  the  injected  flow  and  longitudinal  sound  is  stabilizing  in  the  sense 
that  acoustic  energy  is  absorbed  across  the  injection  region.  Unfortunately,  the  one¬ 
dimensional  model  is  not  capable  of  providing  insight  into  the  details  of  the  absorption 
process.  It  is  of  interest  to  point  out  that  the  "flow-turning"  losses  are  almost  equal  to  the 
increase  in  acoustic  energy  due  to  the  mean  flow.  That  is,  the  sum  of  the  "flow-turning” 
losses  and  the  convection  gains  almost  cancel  with  the  resulting  absorption  of  order  Minj2. 
It  is  important  to  understand  that  in  rocket  motor  stability  calculations,  the  convective  gain 
in  acoustic  energy  really  represents  an  acoustic  energy  loss  because  the  acoustic  energy  is 
convected  out  of  the  the  rocket  interior.  The  important  point  to  be  made  here  is  that  the 
"flow-turning"  losses  are  of  the  same  order-of  -magnitude  as  the  convective  losses. 

Equation.  (57)  also  shows  that  the  "flow-turning"  sound  absorption  is  independent 
of  the  sound  frequency.  This  is  what  one  would  expect  in  a  one-dimensiona!  rigid  wall 
analysis.  These  findings  are  in  agreement  with  the  500  Hz  and  1  KHz  data  presented  in 
Figures  23  &  24.  The  effects  of  both  injection  Mach  number  and  sound  frequency  were 
found  to  be  negligible.  This  behavior  is  not  true,  however,  for  the  2  KHz  data.  Here,  the 
amount  of  sound  energy  that  is  absorbed  increases  almost  linearly  with  M,nj  providing  Mmj 
>0.1.  The  data  suggests  that  the  sound  absorption  behavior  may  be  controlled  by  the  local 
wall  admittance.  Referring  to  Figure  27,  it  may  be  no  coincidence  that  the  values  of  the 
real  part  of  the  wall  admittance  are  almost  equal  at  500  Hz  (AR  =  0.54)  and  1  KHz  (AR  = 
0.59)  while  the  corresponding  value  at  2  KHz  is  much  higher  (AR  =  0.77). 

Now  consider  what  happens  for  finite  admittance  side  walls.  Referring  to  Eq.  (54), 
the  one-dimensional  model  predicts  that  the  interaction  between  transversely  injected  flow 
and  longitudinal  sound  results  in  a  net  gain  of  sound  energy.  Inserting  the  values  of  Ar 
displayed  in  Figure  27  into  Eq.  (54),  the  predicted  predicted  values  are  compared  to 
measured  data  in  Figure  28.  The  predicted  gains  of  energy  are  sufficiently  small  that  the 
effects  of  frequency  variation  can  not  be  seen  -  they  all  collapse  into  the  single  curve 
siiown.  This  slight  increase  is  caused  by  the  finite  values  of  Ar.  It  is  clear  that  the  one- 
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dimensional  model  summarized  by  Eq.  (54)  does  not  agree  with  the  data  displayed  in 
Figure  28. 


Injection  Mach  Number  -  M(inj) 


Figure  28.  Comparison  Between  Predicted  and  Measured  AEE 


The  reader  is  reminded  that  Eq.  (54)  is  valid  only  for  travelling  waves.  However, 
this  restriction  is  not  particularly  severe  because  most  of  the  sound  energy  is  propagating 
downstream.  This  is  evident  from  the  measured  small  values  of  reflection  coefficients  R 
displayed  in  Table  2.  Recall  that  R2  is  a  measure  of  the  relative  fraction  of  the  sound 
energy  propagating  upstream. 


Transient  Solution 

Equations  (35  &  36)  represent  the  wave  equations  governing  the  behavior  of  the 
sound  field  within  the  injection  test  section.  Transient  solutions  to  these  equations  are 
sought  in  order  to  compare  them  directly  to  Culick’s  solution  defined  by  Eq.  (1).  Starting 
with  Eq.  (35),  the  solution  to  the  lowest  order  or  unperturbed  (i.e.,  =  0)  wave 

equation  is  sought  having  the  form 


p0(x  ,t)  =  Cfi' 


<n  (1  ♦  ia)t  -  kx 


=  C(p-a,e‘ ,a*  1 


A  (I  ) 


(58) 
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Equation  (58)  has  been  written  to  show  explicitly  how  the  amplitudes  A(t)  of  the  travelling 
waves  decay  in  time  for  a  >  0.  Substituting  Eq.  (58)  into  Eq.  (35)  yields  the  following 
approximate  solution 


P0(*  -f  )=C0e 


+  i  [cot  —  kx  ] 


(59) 


The  decay  constant  a  is  proportional  to  the  real  part  of  the  wall  admittance  Ar.  This  is 
reasonable  because  the  transient  solution  defined  by  Eq.  (59)  follows  immediately  from  the 
acoustic  steady-state  solution  defined  by  Eq.  (37)  provided  x  and  t  are  connected  using  the 
phase  speed  c  of  the  travelling  wave  as  shown  below 


c  =  —  ->  x  =  ct  (60) 

dt 

The  solution  defined  by  Eqs.  (59  &  60)  assumes  that  the  loss  of  sound  energy  at 
the  wall  is  sufficiently  weak  that  the  sound  wave  can  still  be  characterized  as  plane-wave. 
This  is  consistent  with  small  values  of  Ar.  A  simple  connection  exists  between  the  steady- 
state  solution  defined  by  Eq.  (37)  and  the  transient  solution  defined  by  Eq.  (59).  Equation 
(37)  estimates  the  spatial  decay  of  the  amplitude  of  the  sound  wave  at  a  fixed  moment  in 
time,  while  Eq.  (59)  estimates  the  time  dependent  decay  of  the  amplitude  of  the  sound  wave 
at  a  fixed  position  in  space.  Assuming  representative  cold  flow  test  values  for  the 
parameters  governing  a  of  c  =  34,400  cm/sec,  H  =  6.35  cm  and  Ar  =  0.07  yields  a  = 
379/sec1.  Finally,  if  the  wall  is  rigid,  then  Ar  =  0.  and  no  energy  is  absorbed.  This 
recovers  the  idealized  case  of  plane-wave  sound  propagation  in  an  infinite  constant  area 
duct. 

The  next  step  is  to  apply  the  above  ideas  to  seek  transient  solutions  to  Eq.  (36),  the 
wave  equation  perturbed  by  the  injected  flow.  According  to  Culick's  model,  even  if  the 
walls  are  rigid,  sound  energy  should  be  absorbed  by  the  turning  of  the  injected  flow.  To 
explore  this,  Aw  =  0  on  the  RHS  of  Eq.  (36),  leading  to  the  simplified  wave  equation, 


1 


2  2 

c  dt 


dx 


o 

dU 

C-  \ 

H 

: 

_  F  -  T 

pc 


dX 


dxdt 


(61) 


2 k_ 
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O'  +  2kx  ) 
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In  deriving  Eq.  (61),  the  unperturbed  acoustic  velocity  u0  was  connected  to  the  acoustic 
pressure  pu  via  the  acoustic  momentum  conservation  equation. 


A  solution  to  Eq.  (61)  is  sought  of  the  form 


p,(x  ,  t  )  =  [f  (t  )  +  xg  (t  )]C(p' (“*  *  (62) 

Substituting  Eq.  (62)  into  Eq.  (61),  the  real  part  of  the  solution  is 

pfx  +  Jfp0e!  ^-kx)  (63) 


The  composite  solution,  to  order  Mjnj,  follows  by  inserting  the  solutions  defined  by  Eqs. 
(59  &  63)  into  Eq.  (27)  to  yield 


P(x  ,t) 


1 


i  (cot  -  kx  ) 


1 


+  M. 


l]v 


(tat  -  kx  ) 


(64) 


In  deriving  Eq.  (64),  k  was  replaced  by  0)/c,  and  cMjnj  by  vb  where  vb  is  Culick’s 
expression  for  the  rate  of  fluid  injection  at  the  surface. 


A  discussion  of  the  meaning  of  Eq.  (64)  will  be  deferred  until  Culick's  "flow¬ 
turning"  estimate  defined  by  Eq.  (1)  is  reviewed.  Referring  to  the  derivation  presented  in 
the  Appendix,  p0  represents  the  unperturbed  pressure.  For  travelling  wave  sound 
propagation  in  rigid-walled  constant  area  ducts,  p0  is  defined  as 


it  wt  -  k  x  1 


(65) 


where  kt  =  co/c  using  Culick's  notation.  Setting  the  constant  Co  =  1  without  loss  of 
generality,  and  substituting  Eq.  (65)  into  Eq.  (1)  yields  the  following  estimate  of  ctpj- 
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(66) 


'FT 


H 


The  "flow-turning"  decay  constant  defined  by  Eq.  (66)  represents  a  loss  of  energy 
characterized  by  the  expression 


p(x,  t ) 


V  t 


0 

77 


+  /  ( ot  -  kx) 


(67) 


The  acoustic  energy  loss  in  Eq.  (67)  is  exponential.  Now  compare  Eq.  (67)  to  the 
corresponding  energy  loss  expression  defined  by  Eq.  (64).  If  one  assumes  that  Minj  «  1, 
then  the  term 


VA 

H 


M  TTn 

*»  H  JJ 


=  e 


■“»w) 


(68) 


Using  Eq.  (68),  Eq.  (64)  may  be  approximated  by  the  expression 


p(x  ,t)  = 


,  tM  X  )  -  kx  ) 

+  Mti  ^')+  U0"  ) 


(69) 


A  comparison  between  Eq.  (67)  and  Eq.  (69)  verifies  Culick’s  model  providing  the 
approximation  defined  by  Eq.  (68)  is  valid.  This  is  important.  The  implication  here  is  that 
if  enough  terms  were  solved  for  in  the  regular  perturbation  solution  defined  by  Eq.  (27), 
then  their  sum  would  equate  to  the  exponential  expression  derived  by  Culick.  It  is  clear 
that  this  is  a  time  consuming  and  laborious  task.  This  is  a  dramatic  way  of  illustrating  the 
cleverness  of  Culick's  approach  wherein  he  sought  solutions  only  for  the  perturbed 
wavenumbers  and  not  for  the  perturbed  pressures;  a  far  easier  task. 
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It  is  of  interest  to  compare  the  relative  importance  between  the  "flow-turning"  time 
decay  constant  defined  by  Eq.  (69)  and  the  finite  admittance  wall  absorption  time  decay 
constant  defined  by  Eq.  (59).  The  ratio  of  the  decay  constants  are 


TT 

'  V  '  4,  <70> 

H 


Substitution  of  representative  test  values,  Mjnj  =  0.003  and  AR  =  0.07,  into  Eq.  (70)  yields 


=  0.04 


(71) 


Equation  (71)  shows  that  even  at  the  highest  injection  flow  rate  tested,  the  sound 
energy  losses  due  to  "flow-turning"  are  much  smaller  than  the  sound  energy  absorbed  by 
the  finite  admittance  side  walls.  Using  the  above  parameters  as  an  example,  otpr  =  16  sec'1 
and  aAw  =  379  sec1.  Both  the  data  summarized  in  Figures  23  &  24  and  the  above  analysis 
show  that  the  "flow-turning"  losses  are  much  smaller  than  the  losses  due  to  the  finite 
admittance  side  walls. 


CONCLUDING  REMARKS 


Test  results  show  that  acoustic  energy  is  absorbed  by  the  transverse  injection  of 
ambient  air  into  an  acoustically  energized  duct.  The  data  can  be  partitioned  into  low  and 
high  frequency  behaviors.  For  frequencies  below  1  KHz,  the  energy  absorption  was 
negligible  and  independent  of  frequency  and  injection  flow  rate.  At  2  KHz,  just  below  the 
first  transverse  cut-on  mode  of  the  duct,  over  1  dB  of  sound  was  absorbed  -  a  significant 
amount  -  at  an  injection  Mach  number  rate  of  0.003.  In  all  test  cases,  the  sound  absorbed 
by  the  finite  admittance  side-walls  was  significantly  larger  than  that  absorbed  by  the 
injected  flow. 

The  effects  of  refraction  of  the  sound  pressure  by  the  mean  flow  velocity  gradients 
were  found  to  be  negligible.  This  was  not  surprising  in  light  of  the  very  small  injection 
flow  rates  used  in  these  tests. 
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One-dimensional  steady-state  and  transient  "flow-turning”  solutions  were  derived. 
For  the  special  case  of  rigid  walls  and  low  injection  rates,  the  steady-state  solutions  showed 
that  the  "flow-turning"  losses  were  approximately  equal  to  10Log(l-Mav).  These  losses 
were  almost  cancelled  by  the  convected  energy  gains  which  were  equal  to  10Log(l+Mav). 
The  contribution  of  the  Finite  admittance  side  walls  resulted  in  overall  very  small  acoustic 
energy  gains.  The  one-dimensional  solutions  inadequately  predicted  the  measured  acoustic 
energy  losses. 

The  transient  solution  compared  very  favorably  with  Culick's  model.  Both 
estimated  almost  the  same  decay  constants.  Using  the  data  as  input,  the  models  showed 
that  the  decay  constants  associated  with  the  absorption  of  sound  by  the  injected  flow  were 
much  smaller  than  the  decay  constants  associated  with  the  absorption  of  sound  by  the  Finite 
admittance  side  walls.  These  Findings  cannot  be  blindly  used  to  infer  that  "flow-tuming" 
losses  are  unimportant  in  rocket  motor  applications.  They  are  valid  only  for  the  cold  flow 
test  conditions  of  the  facility.  For  example,  the  values  of  the  wall  admittance  within  rocket 
motors  are  generally  much  smaller  than  the  admittance  of  the  Glasstrate  material  used  in  this 
test  program.  An  order-of-magnitude  estimate  of  typical  solid  rocket  admittance  values  can 
be  determined  from  response  function  measurements.  The  response  function  and  the  wall 
admittance  are  connected  as  shown  below, 


m' 


P 

According  to  Dr.  R.  S.  Brown12,  a  typical  value  of  R  is  about  2.  Combining  this  with 
Mjnj  =  0  .003  leads  to  Aw  =  0.006  which  is  an  order-of-magnitude  less  than  the  Glasstrate 
admittance  values.  Substituting  these  values  into  Eq.  (70)  yields  the  following  estimate  of 
the  importance  of  Culick's  "flow-tuming"  losses  relative  to  the  losses  generated  by  liner 
walls  within  solid  rocket  motors  , 


M 


-  1 


(72) 


a 


FT 


a 


=  1.  5  -  0  (1) 


(73) 


It  was  shown  in  Section  4.2  that  the  one-dimensional  model  was  inadequate  in 
predicting  the  measured  acoustic  energy  losses  across  the  injection  test  section.  Referring 
to  Eqs.  (66  or  69),  the  one-dimensional  model  shows  the  "flow-tuming"  losses  to  be 
inversely  proportional  to  the  duct  height  H.  This  makes  sense  because  H  is  the  only 
relevant  length  scale  available.  Two  (or  three)  -  dimensional  models  may  show  that  the 
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relevant  length  scale  over  which  most  of  the  "flow-turning"  losses  occur  may  be 
significantly  less  than  H.  The  idea  here  is  that  near  the  head  end,  the  flow  penetrates 
efficiently  into  the  rocket  center.  Near  the  exhaust  end,  however,  the  penetration  is  less 
efficient.  This  suggests  an  overall  effective  length  scale  less  than  H.  If  this  proves  to  be 
true,  then  H  in  Eqs.  (66  or  69)  would  be  replaced  with  a  smaller  term  resulting  in  an 
enhancement  of  the  contribution  of  the  "flow-turning"  losses  to  the  stability  of  solid  rocket 
motors.  The  poor  agreement  between  predicted  and  measured  acoustic  energy  losses 
underscores  the  need  to  develop  more  accurate  models  of  the  acoustic  energy  losses 
associated  with  the  injection  of  combustion  products  into  the  interior  of  solid  rocket 
motors. 


One  of  the  important  findings  of  this  investigation  is  that  it  clearly  illustrates  the 
difficulty  of  using  cold  flow  laboratory  facilities  to  measure  "flow-turning"  losses.  Most  of 
the  sound  absorption  takes  place  at  the  finite  admittance  side  walls  through  which  the  mean 
flow  is  injected.  Experimentally,  this  translates  into  trying  to  accurately  measure  small 
"flow'-tuming"  losses  as  the  differences  between  two  relatively  large  terms.  One  needs  to 
think  of  clever  ways  to  get  around  this  problem. 


52 


REFERENCES 


1.  Culick,  F.E.C.,  "The  Stability  of  One-Dimensional  Motions  in  a  Rocket  Motor”, 
Combustion  Science,  and  Technology.  Vol.  7,  pp.  93-97,  1977. 

2.  Culick,  F.E.C.,  "Stability  of  Three-Dimensional  Motions  in  a  Combustion  Chamber", 
Combustion  Science,  and  Technology.  Vol.  10,  pp.  109,  1975. 

t  3.  Culick,  F.E.C.,  "Remarks  on  Entropy  Production  in  the  One-Dimensional 

Approximation  to  Unsteady  Flow  in  Combustion  Chambers",  Combustion  Science,  and 
Technology,  Vol.  15,  pp.  93,  1977. 

4.  Culick,  F.E.C.,  Magiawala,  K.,  Wat,  J.,  Awad,  E.,  and  Kubota, T.,  Measurements  of 
Energy  Losses  Associated  with  Interactions  between  Acoustic  Waves  and  Non-Uniform 
Steady  Row.  AFRPL-TR-81-22.  July.  1981. 

5.  Flandro,  G.A.,  "Solid  Propellant  Acoustic  Admittance  Corrections",  Journal  of  Sound 
and  Vibration.  Vol.  36,  No.  3,  1974,  pp.  297-312. 

6.  Vuillot,  F.  and  Kuentzmann,  P.,  "Flow  Turning  and  Admittance  Correction:  An 
Experimental  Comparison",  Journal  of  Propulsion  and  Power.  Vol.  2,  No.  4,  July-August 
1986,  pp.  345-353. 

7.  Baum,  J.D.  and  Levine,  J.N.,  Acoustic-Mean  Flow  Interaction.  AFRPL  TR-86-065, 
January,  1987. 

8.  Van  Moorhem,  W.  K.  ."Theoretical  Basis  of  the  Flow  Turning  Effect  in  Combustion 
Stability",  UTEC  ME  80-033,  Univ.  of  Utah,  Salt  Lake  City,  Utah.,  March,  1980. 

9.  Culick,  F.  E.  C.,  "Rotational  Axisymmetric  Mean  Flow  and  Damping  of  Acoustic 
Waves  in  a  Solid  propellant  Rocket",  AIAA  Journal.  Vol  4,  No.  8,  August,  1966,  pp. 
1462-1464. 

10.  Hersh,  A.  S.  and  Walker,  B.,  Experimental  Investigation  of  Rocket  Motor  Flow- 
Turning  Acoustic  Losses.  AFRPL  TR-84-009,  May,  1984. 

11.  Pridmore-Brown,  D.C.,  "Sound  Propagation  in  a  Fluid  Flowing  Through  an 
Attenuating  Duct",  J.  Fluid  Mech..  Vol. 4,  pp.  1393,  1958. 

t  12.  Brown,  R.  S.,  Telephone  communication  concerning  solid  rocket  response  functions, 

Sr.  Scientist,  United  Technologies  Chemical  Systems,  San  Jose,  CA.,  12  January,  1988. 


53 


APPENDIX  -  DERIVATION  OF  FLOW  TURNING  LOSSES 


For  the  convenience  of  the  reader,  a  detailed  derivation  of  Eq.  (1)  is  derived  below 
in  order  to  understand  what  is  meant  by  the  term  "flow-turning"  losses.  The  derivation  is 
based  in  part  upon  the  earlier  work  of  Culick,  et.  aH.  Consider  the  control  volume 
sketched  below  consisting  of  the  transverse  injection  of  flow  into  a  uniform  channel  of 
height  H  with  no  combustion, 


Injected  Mass  Flow  Rate 


Injected  Mass  Flow  Rate 


Figure  A-1.  Control  Volume  Used  in  Model  Derivation 


The  injected  flow  is  assumed  to  be  inviscid.  Restricting  the  flow  within  the  channel 
to  be  one-dimensional,  the  conservation  equations  of  mass  and  momentum  flux  can  be 
written 


dp 

at 
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J<-(pU)=  gpv 


and 
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The  derivation  of  the  above  equations  assumes  that  mass  is  transversely  pumped 
into  the  control  volume  at  the  rate  pv  per  unit  area.  The  usual  practice  is  assumed  that  the 
instantaneous  fluid  properties  can  be  decomposed  into  mean  and  acoustic  components, 

p  =  p+  p’;  p  =  p+  p’;  u  =  U(x)+  u';  v  =  V  +  v'  (A-3) 


The  mean  velocities  U  &  V  are  connected  via  the  steady-state  conservation  of  mass  flux 
equation.  Inserting  Eq.  (A-3)  into  Eqs.  (A- 1  &  A-2)  yields  the  following  conservation  of 
acoustic  mass  and  momentum  flux  equations, 


§r +  plx =  H(Vp  +  py,)  ■  e + wp' 


and 


djr^dpi  2r 

P  dt  dx  H 


pUv ’  +  V  fpu  '  +  Up' 


a,du'  dU 

-pvw-p*r 


(A-4) 


(A-5) 


The  RHS  of  Eqs.  (A-4  &  A-5)  contain  five  terms  which  were  generated  from  the 
injected  flow.  Culick  observes  that  the  RHS  terms  may  be  considered  to  be  source  terms 
for  the  LHS  by  expanding  the  acoustic  pressure,  density  and  particle  velocity  in  powers  of 
V,  the  mean  injection  speed.  With  regard  to  deriving  the  "flow-turning"  losses  identified  as 
Eq.  (1)  in  the  Introduction,  the  only  (source)  term  to  be  retained  is  the  one  identified  within 
the  brackets  on  the  RHS  of  Eq.  (A-5).  The  resulting  simplified  equations  are 


dp’  _  du' 

ar  +  pa^r 


=  0 


(A-6) 


for  the  continuity  flux  equation  and 
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for  the  momentum  flux  equation.  Observe  from  comparing  Eqs.  (A-6  &  A-7)  with  Eqs. 
(A-4  &.  A-5)  that  the  "flow-turning"  losses  represent  only  one  of  many  acoustic  sources. 
The  relative  importance  of  this  source  is  compared  to  the  others  in  Section  4. 

Since  no  combustion  is  assumed  to  be  taking  place  and  the  rate  at  which  air  is 
injected  into  the  channel  is  very  low,  fluctuations  of  the  fluid  pressure  and  density  take 
place  adiabatically  and  the  acoustic  energy  flux  conservation  equation  may  be  written 


p'  =  cy 


(A-8) 


where  c  is  the  fluid  speed  of  sound.  The  wave  equation  is  derived  in  the  usual  way  by 
using  Eq.  (A-8)  to  eliminate  p'  and  differentiating  Eq.  (A-6)  with  respect  to  time  and  Eq. 
(A-7)  with  respect  to  distance  to  yield 


^P!  _  _  _  IgySul 

ax’  c!  at1  hp  ax  <a-9> 

In  the  absence  of  injected  flow,  Eq.  (A-9)  reduces  to  the  familiar  wave  equation, 


aV.  i  aV, 

ax2  c2  at2 


(A- 10) 


where  p#  denotes  that  the  sound  pressure  is  unperturbed  by  the  injected  flow. 

The  boundary  conditions  associated  with  Eqs.  (A-9  &  A- 10)  require  special 
comment.  Culick  assumes  the  constant  area  channel  is  bounded  by  end  walls  at  x  =  0,L. 
The  unperturbed  acoustic  particle  velocities  at  these  locations  vanish.  Applying  these 
boundary  conditions  to  the  unperturbed  and  perturbed  wave  equations,  yields,  respectively 


=  0 


(*  =0,  L) 


and 


(All) 
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(x  =0,  L) 


(A-12) 


cp'  2  pVu\ 

A"  "  H 


The  boundary  conditions  specified  by  Eq.  (A-12)  are  not  clear.  We  presume  that  Culick 
allows  the  small  amount  of  volume  flow  generated  within  the  channel  to  "somehow"  leak 
through  the  end  walls  at  x  =  0  and  x  =L;  otherwise  u'  vanishes. 

The  wave  equations  defined  by  Eqs.  (A-9  &  A- 10)  are  simplified  by  assuming  the 
standard  time  dependence 


p'(x,  t)  =  ef“‘p(x);  u'(x,t)  =  e'“'u(x) 


(A- 1 3) 


Substituting  Eq.  (A- 13)  into  Eqs.  (A-9  &  A- 10)  results  in  the  following  simplified 
expressions. 


dx 


Hr  dx 


(A- 14) 


and 
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(A- 15) 


Culick  used  a  rather  clever  scheme  to  derive  Eq.  (1)  of  the  Introduction.  To  begin 
with,  he  multiplied  Eq.  (A-14)  by  p0  and  Eq.  (A-13)  by  p,  subtracted  one  from  the  other 
and  integrated  over  the  length  L  of  the  channel  to  yield 


J 


ft 


(A-16) 


Integrating  the  first  by  parts  and  applying  the  boundary  conditions  defined  by  Eqs.  (A-l  1 
&  A-12)  produces  the  following  expression  fork2. 
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where 

L  L 

E,  =  jpp.dx  «  Jp2#dx 

•  • 


(A- 17) 


(A- 18) 


Now  Eqs  (A-7  &  A- 13)  connect  the  velocity  u  to  the  pressure  p.  To  lowest-order, 


iptou  =  - 


dp 

dx 


(A-19) 


Substituting  Eq.  (A-19)  into  Eq.  (A- 18)  yields 
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The  LHS  of  Eq.  (A- 17)  can  be  approximated  by  the  expression 


k  -  k2  =  (  — — -  k  s  k  -  2ik  —  -  k*  =  -  2 ik  — 
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(A -21) 


The  derivation  closes  by  combining  Eqs.  (A-21  &  A-20)  yielding  the  desired  result 
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(A-22) 


Referring  to  (Eq.  1)  in  the  Introduction,  the  quantity  q  represents  the  perimeter  over 
which  fluid  is  injected  into  the  channel  interior  of  cross-section  Sc.  Referring  to  Figure  A- 
1,  q  =  2W  where  W  is  the  width  of  the  channel.  The  channel  cross-section  area  Sc  =  WH. 
Thus  q/Sc  =  2/H.  Finally,  the  quantity  mb  =  pV.  With  these  understandings,  Eqs  (1  &  A- 
22)  are  equivalent 
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